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Abstract 

Using the Embedded Atom Method as developed by Voter and Chen in combination with the variable 
metric /quasi- Newton and our own Aufbau/Abbau methods, we have identified the three most stable 
isomers of Auat clusters with N up to 150. For the first time clusters with tetrahedral symmetry are 
found to form the ground states of Auir and AU34. The AU54 icosahedron without a central atom and 
the Aui46 decahedron are found to be particularly stable, whereas the highly symmetric second and third 
Mackay icosahedra that could have been obtained for A^ = 55 and 147, respectively, do not correspond to 
the particularly stable structures. The three lowest-lying isomers of AU55 and AU147 are low-symmetrical 
structures. Various structural and energetic properties are analysed, such as stability function, occurrence 
of magic-sized clusters, construction of icosahedral and fee shells, and cluster growth. 
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This work is dedicated to Prof. Dr. Wolf Weyrich, Physical Chemistry, University of Konstanz, Germany, 
on the occasion of his 65th birthday. 



1 Introduction 

Since the first part of tlie 20tli century it has been recognized that theoret- 
ical studies can constitute an important ingredient of science, first through 
the dedication of chairs to the field of theoretical physics and later through 
the establishment of chairs in theoretical chemistry. Simultaneously, by 
defining 'theory' as an independent part of science, the interactions be- 
tween theory and experiment have in many cases been reduced. Persons 
who unite these two parts are, therefore, of immense importance and, si- 
multaneously, able to contribute to science in a way that is unmatched by 
most colleagues. A person, occupying a chair of physical chemistry, and 
being active in both experimental and theoretical studies of the properties 
of matter, is, accordingly, unique. Such a person is Wolf Weyrich, who 
has made significant contributions to the understanding of the properties 
of a large range of systems, going all the way from smaller molecules to 
extended solids. In the later years his interests have also turned towards 
the properties of nanoparticles that somehow lie in between the materials 
of his earlier interests. Metal nanoparticles are thereby of central interest. 
Therefore, it is with pleasure that we dedicate the present work to Wolf 
Weyrich on the occasion of his 65th birthday. 

Clusters form an important link between isolated atoms and molecules 
at one extreme and bulk solids at the other. Their large surface-to- volume 
ratio gives them unique physical and chemical properties. On the other 
hand, the combination of finite with large size makes it difficult to char- 
acterize and analyze their properties in detail. From a theory point of 
view the central difficulty lies in the determination of the structure of 
the system of interest and only through comparison with experimental 
information more definite statements about their properties can be made. 

Gold clusters represent some of the mostly studied clusters (for a 
detailed discussion, see Ref. E]). They have recently been investi- 
gated in connection with the synthesis of nanostructured materials and 
devices |21 HI El El C]- Their structural and energetic properties have 
been studied with High-Resolution Electron Microscopy (HREM) and var- 
ious spectroscopic techniques 0131101111111211131111211111111111111 
1201 121 22, 23 24 25 . Literature concerning small Auat clusters is en- 
riched with numerous investigations based on density-functional methods 

ElElElOlElEIlElElElElElEZlElEllSllIIlEllllllll that 

are not yet capable of giving a definite answer to the problem at what 
cluster size the structural 2D - 3D transition occurs. Recent studies com- 
bining theory and experiment |31 01 E| |3SI show that the gold clusters 
are planar at least up to A^ = 7 for Ref. 0H|, or up to A^ = 12, according 
to Hakkinen 45 and Furche et al. |47j . 

However, global structure optimization is difficult when using ah initio 
methods already at very small cluster sizes. Nevertheless, some studies 
in this direction exist. Thus, the authors of Ref. j^E^I performed den- 
sity functional calculations on clusters containing more than 30 atoms. 



relaxing selected high-symmetric configurations. Alternatively, the global 
optimizations of larger clusters are all based on approximate methods 
like molecular dynamics IllEniEIlE21E3ElESlEniE3EHlEniEni 
and semiempirical potentials like the EAM |13[ I61L I62| . Sutton-Chen 
[S5| , Murrell-Mottram |^ ESI ; or the many-body Gupta potential |^ 
IS7I EHl EHl EOl [^' Using these methods, unbiased structure optimiza- 
tions were performed up to the 80-atom cluster. Medium-sized clusters 
( 80 < TV < 150) have hardly been studied. Besides the first-principles 
study of Haberlen et al. j^T] and the EAM calculations by Cleveland 
et al. ^1 1^ ins considering particular structural motifs, there exists 
essentially no further investigation on the clusters in this size range. 

In most of the studies, special attention is paid to the so-called 'magic- 
numbered' clusters, that possess closed electronic and/or geometric shells. 
Various studies on the smallest 'magic' cluster Auia have identified the 
formation of an icosahedron [211 123 EHl 1^ IMl IS] ■ Only the authors of 
Ref. |34[I35) found a disordered structure as the lowest-lying isomer for this 
cluster size. On the other hand, semiempirical potentials ,6^ ESI |^ [^ 
and the density functional study by Hakkinen et al. j^H] on the Auas 
cluster predict the truncated octahedron to be the global minimum for this 
cluster size. However, on the basis of first-principles and Gupta potential 
calculations, the authors of Ref. |HEI EOI state that a disordered structure 
is actually lower in energy than the symmetric. Ultimately, it may be 
suggested that the obtained structure depends sensitively on the type of 
the potential, since by using another form of the same potential, Darby et 
al. j71j found a truncated octahedral structure to be the global minimum 
of Au38. The situation is more clear for the Auss and AU75 clusters, where 
a disordered structure [031 1001 lOZl JOHl [ZOl [^ or a Marks decahedron (m- 
Ds/i) |67[ 165] seem to constitute the global minima. 

Several groups have performed calculations on larger clusters by mini- 
mizing the total energy of initially chosen symmetric structures, although 
it may be feared that the structures may not be those of the global total- 
energy minima. Moreover, only few cluster sizes were studied — the 
octahedral Auyg j^ 103 and Aui4o IOHj the decahedral Auioi, Aung, 
and Aui46 |02> and the icosahedral AU147 |2Z1|^. The structures and 
energetics of the clusters between these high symmetrical ones remain 
scarcely investigated. 

The purpose of the present study is to carry out unbiased calculations 
on small and intermediate gold clusters, and, subsequently, to investigate 
the occurrence of magic clusters and the growth patterns, especially for 
larger clusters. To our knowledge, two earlier studies on gold clusters 
with the same version of EAM have been performed previously, i.e., the 
work of Rey and coworkers |72[ [75] and that of Sebetci and Giiveng [TJ . 
The first study considered only the energetics and stability of small gold 
clusters (2 < A^ < 23). Sebetci et al. used a basin- hopping Monte Carlo 
minimization approach to find the global minima of AIat, Au^r, and Pt^r 
clusters with N < 80. The total energies, point groups, and structural 



assignments were presented. 

In the present study the structure and energetics of the three most 
stable isomers of small and medium-sized Auat clusters with 2 < iV < 150 
have been determined for each cluster size by using a combination of the 
embedded-atom method in the version of Voter and Chen J78L 1791 150] , 
the variable metric /quasi- Newton method, and our own Aufbau/Abbau 
method. The paper is organized as follows. In Sec. 2 we briefly outline 
the embedded-atom method, and in Sec. 3 we present our structural- 
determination methods. The main results are given in Sec. 4, and a brief 
summary is offered in Sec. 5. 

2 The Embedded Atom Method 

The main idea of the EAM was initially presented by Daw, Baskes, and 
Foiles (DBF) [73|7i[77| in 1983-1986, and since then the generality of 
the functions of the EAM of DBF has been successfully tested through 
numerous applications to different systems of metals and alloys, including 
defects, surface and interface structures, surface and bulk phonons, etc. 
In a previous paper |HJ| we reported results for the global minima of 
Nijv, CuAT, and Au^v clusters with up to 60 atoms, obtained with two 
different versions of the EAM. There we discussed the incapability of the 
so-called DBF version of EAM to describe properly the properties of the 
smallest gold clusters, which could be related to the parameterization of 
the potential only to bulk properties. The version developed by Voter 
and Chen fTSIEHUHOI takes into account also the properties of the dimer, 
which makes this method more suitable for the description of the smallest 
clusters. Accordingly, in this study we use this version of the EAM for 
the calculation of the total energy of a given cluster. 

The principle of the method is to split the total energy of the system 
into a sum over atomic energies: 

N 
i 

The embedding energy is obtained by considering each atom as an im- 
purity embedded into a host provided by the rest of the atoms. The 
electron-electron interaction is presented as sums of short-ranged, pair 
potentials. Accordingly, 

i i^j 

where p'' is the local electron density at site i, Fi is the embedding energy, 
i.e., the energy required to embed an atom into this density, and 4>ij is 
a short-range potential between atoms i and j separated by distance r^ . 



The pair potential, according to the Voter-Chen version, is taken to be a 
Morse potential, 

0(r) = Dm[1 ~ e-«M('--i?.M)]2 _ ^^^ (3) 

where the three parameters. Dm, Rm, and um, define depth, position of 
the minimum, and a measure of the curvature at the minimum, respec- 
tively. The local density at site i is assumed being a superposition of 
atomic electron densities: 

where p" (r^ ) is the spherically averaged atomic electron density provided 
by atom j at the distance r^ . The density function is taken as the density 
of a hydrogenic 4s orbital: 

p(r)=r6[e-'^'- + 29e-2/5n (5) 

where /? is an adjustable parameter. Because r^e'^"^ turns over at short r, 
the second term has been added to maintain the monotonically decreasing 
character of p{r) at shorter r. This 4s orbital density, appropriate for Ni 
and Cu, also works well for gold. To ensure that the interatomic potential 
and its fist derivatives are continuous, both 0(r) and p{r) are cut off at 
'i'='>'cut- In the fitting procedure, the five parameters defining (f)(r) and 
p{r) {Dm, Rm, ocm, 13 m and rcut) are optimized by minimizing the root- 
mean-square deviation between the calculated and reference properties 
of different selected systems like molecules, surfaces, solids, and defects. 
Because F{p^) is redefined for each choice of the parameters, the potential 
always gives perfect agreement with experimental values such as oq, i?coh 
and the bulk modulus B. The reference properties are the three cubic 
elastic constants (Cu, C12 and C44), the unrelaxed vacancy formation 
energy {El^^), and the bond length {Re) and bond strength {D^) of the 
diatomic molecule. The values of pf, (j)ij and Fi{pi) that were used by the 
Voter-Chen version, are available in numerical form for Ni, Pd, Pt, Cu, 
Ag, Au and Al. 

Our reason for choosing the EAM was dictated by the good agreement 
to experiment, as well as to first principles calculations, and last but 
not least by the high computational efficiency allowing one to investigate 
clusters with more than 100 atoms without severe constraints on the initial 
geometry, which is impossible with first principles methods. In a previous 
work 81 we performed calculations on smaller gold clusters with 2 < N < 
60 atoms comparing the EAM of Daw, Baskes, and Foiles (DBF) and the 
version used in this work, and we found that the DBF overestimated 
the binding energy of the dimer by 209% and underestimated the bond 
distance by 37%. For comparison, the EAM of Voter and Chen gives 
dimer binding energy corresponding to 99.6% of the experimental value. 



and bond distance that is 92.2% of the experimental value. For this reason 
we chose to work with the Voter-Chen version that describes correctly the 
dimer properties. 

3 Structure optimization 

Using expression |^ we can calculate the total energy of any cluster with 
any structure as a function of structure, i.e., of the atomic coordinates 
{Ri}, Etot{Ri,R2, ■■■ ,Rn)- In order to obtain the closest local total- 
energy minimum we use the variable metric/quasi-Newton method 1881 . 

For searching the global minima we have developed our own Auf- 
bau/Ahhau method that is described in details in previous works |MII55] . 
It consists of the following steps: 

1) We consider two cluster sizes with N and N + K atoms with K ~ 
5 — 10. For each of those we randomly generate and relax a large set of 
structures, from which those with the lowest total energy are selected. 

2) One by one, each of the N atoms is displaced randomly, and the 
closest local minima is determined. If the new structure has a lower 
total energy than the original one, this new one is kept, and the old one 
discarded. This is repeated approximately 1000 — 2000 times depending 
on cluster size. 

3) This leaves us with two 'source' clusters, Auat and Aun+k with 
their lowest total energies. One by one an atom is added at a random 
position to the structure with A'' atoms (many hundred times for each 
size), and the structures are relaxed. In parallel, one by one an atom is 
removed from the structure with N + K atoms — for each intermediate 
cluster with TV' atoms we consider all A^' -I- 1 possible configurations, that 
one can obtain by removing one atom from the Aujv'+i cluster. From 
the two series of structures for A^ < M < A^ -I- if those structures of the 
lowest energies are chosen and these are used as seeds for a new set of 
calculations. First, when no lower total energies are found in the two sets 
of calculations, it is assumed that the structures of the global-total-energy 
minima have been identified. 

4 Results and Discussion 

4.1 Small gold clusters 

Like all other semiempirical potentials, the one of the embedded-atom 
method does not include explicitly the electrons and their orbitals. There- 
fore, such a potential tends to prefer high-symmetry, compact structures, 
whereas structures of lower symmetry that can be explained through elec- 
tronic effects are not found. As a result, for the smallest gold clusters, 
where spin-orbit interactions play an important role, our global-minima 
structures are compact (see Table ^, and the planar structures (that are 



believed to be those of the true total-energy minima) only metastable. In 
the size range iV = 4 - 7 first principles studies obtain these 3D configu- 
rations as higher-lying isomers, which can serve as an example of how the 
inclusion of electronic effects can change the energetic ordering of the iso- 
mers. On the other hand, the addition of electronic effects in the semiem- 
pirical potentials would restrict their use only to small and relatively larger 
clusters with pre-chosen structures. An appropriate choice in this respect 
could be the Density Functional Tight-Binding methods (DFTB) that 
include explicitly the electrons and are computationally more efficient 
than the common density functionals. Actually, in a recent study |82| 
we demonstrated the important role that the electronic effects can play 
for the binding energy and the stability functions. However, except for 
A^ = 4, where a rhombus was the lowest-lying isomer according to the 
DFTB method, all the ground state structures for the smallest gold clus- 
ters Au5 - Aug had 3D shapes. On the other hand, even the most recent 
density-functional studies are still not in agreement at which cluster size 
the structural transition 2D - 3D occurs. According to the LDA study 
of Wang et al. [3S] , it is the pentagonal bipyramid that forms the global 
minimum of AU7. Remade and Kryachko 03| suggested that gold clus- 
ters are planar at least up to iV = 9, while Walker |33| predicted that the 
transition occurs at Aun. Using ion mobility measurements and ab-initio 
molecular dynamics Kappes et al. [31103 found that the 3D transitions 
occur at Auj~2 and Au^. The same group studied the adsorption of CO 
on isolated gold cluster cations in the size range N = 1-65. The smallest 
clusters with 4 < A^ < 6 as well as Aug were found to be planar, while 
for Au7 the global minimum was a 3D structure, but not a bipyramid, in 
contrast to the results of Wang et al. |2S1 ■ In a combined experimental 
and theoretical study Hakkinen and coworkers 0B] confirmed the 2D - 
3D transition at Au^jj however, Xiao and Wang 021 suggested that for 
the neutral clusters this transition occurs first at Auis. In niost of the 
cases the planar structures are competing with 3D isomers, and the ener- 
getic differences are insignificant, which in turn means that the ordering 
of the isomers depends strongly on the used functional and the starting 
conditions. For example, Hakkinen et al. |33| compared the global min- 
ima of relativistic and nonrelativistic Au^ clusters and found that for the 
nonrelativistic gold the lowest-lying isomer was a capped octahedron that 
corresponds to our second isomer for this cluster size. At larger cluster 
sizes the potential used in this study yields results in agreement with 
density-functional and experimental studies. The study of Hakkinen et 
al. 1231 on the Ausg cluster predicted the truncated octahedron to be 
the global minimum, and a recent experiment |25| showed that the Au^g 
cluster most probably is not an icosahedron, but a structure with a low 
symmetry, in agreement with our results. 

In summary, we can conclude that although our results for the smallest 
gold clusters correspond to higher-lying isomers within the first principles 
methods, due to the lack of electronic effects, our method is sufficiently 



accurate in describing the larger gold clusters with N > 9, where most 
probably the planar structures begin to compete with 3D configurations. 

4.2 Energetic properties 

The high stability of the so-called 'magic-numbered' clusters has become 
a subject of great interest in connection with its relevance in the medicinal 
and the colloidal chemistry, as well as in the production of catalysts and 
high-tech nanomaterials. 

In Fig. ^ we show the binding energy per atom for the global- minima 
structures, as well as the difference between the total energies of the 
lowest-lying isomers obtained by us and those found by Sebetci et al. |74| 
using exactly the same potential for the interatomic interactions. One can 
see that the latter difference increases almost linearly with the number of 
atoms and has its maximum at iV = 79, where we obtained a truncated oc- 
tahedron in contrast to the structure with D^h symmetry found by Sebetci 
and Giiveng. Except for few cases, the total-energy difference is marginal 
(about 5 meV/atom, which may be due to numerical differences), giving 
support for the quality of both theoretical approaches in optimizing the 
structure. At A^ = 52 they obtained an uncentered icosahedron-like struc- 
ture with C2h symmetry, that lies energetically between our second and 
third lowest isomers. 

In order to identify the particularly stable clusters we have considered 
the following criteria. The clusters can be considered as very stable if their 
binding energy per atom is much larger than that of the two neighboring 
clusters. This can be quantified through the stability function, £'tot(-^ + 
1.1) + EtotiN - 1.1) - 2£;tot(iV.l), where EtotiN.k) is the total energy of 
the energetically fc-lowest isomer of the Auat cluster. This function, that 
has maxima for particularly stable clusters, is shown in Fig. [3 Here we 
can identify a large number of particularly stable clusters, i.e., so-called 
magic clusters. These are found for A^ = 4, 6, 10, 13, 15, 17, 23, 28, 
30, 36, 38, 40, 42, 45, 49, 54, 58, 61, 64, 66, 68, 73, 75, 77, 79, 82, 84, 
89, 92, 95, 101, 109, 111, 116, 118, 124, 128, 133, 135, 140, 144, and 
146. The most pronounced peaks (marked in the figure) occur at A^ = 
13, 30, 40, 54, 75, 79, 82, 124, 133, 140, and 146. In agreement with 
Sebetci and Giiveng [TJ , the 54-atom icosahedron without a central atom 
is found to represent a magic-numbered cluster, whereas the AU55 cluster 
does not. The latter possesses a distorted icosahedral structure with €3^ 
symmetry, lying 0.374 eV lower than the perfect icosahedron, 2.9 eV lower 
than the decahedron, and 3.27 eV lower than the cuboctahedron. In 
our study, all the three lowest-lying isomers of AU55 have lower energy 
than the symmetric structures, in agreement with previous studies where 
disordered configurations were found as global minima for AU55 |66l 1671 
COlIZIl- For Au38 and AU75, a cuboctahedron 29 , 63 , ^ EH El and a 
Marks decahedron |63l 1671 W^ I74j were obtained, in agreement with first- 
principles and semiempirical studies. However, two studies employing 



the many-body Gupta potential identified amorphous structures as those 
of the global minima of Auas HHI IZOli which is most probably due to 
the parameterization of the potential, since Darby et al. |71| found an 
octahedron as the lowest-lying isomer by using another version of the 
same potential. 

Another striking result of our study is that the 146-atom Marks dec- 
ahedron represents a peak in the stability function, whereas the AU147 
icosahedron does not. According to our study, the third Mackay icosahe- 
dron lies 2.89 eV lower than the cuboctahedron, which in turn is 2.53 eV 
lower than the decahedron, but 0.37 eV higher than a disordered structure 
with partly decahedral construction. To our knowledge, this is the first 
study predicting a disordered global minimum for the AU147 cluster. 

According to our other criterion for a particularly stable cluster, such a 
cluster occurs if the energy difference between the two energetically lowest 
isomers EtQt{N.2) — Etot{N.l) is large. This energy difference is shown 
in Fig. |31 and comparing to Fig. [3 we can see that many of the clusters 
that are particularly stable according to the first criterion are stable also 
according to the second one. 

4.3 Structural properties 

In this subsection, instead of discussing in particular the structures of 
the individual clusters, we shall introduce different quantities that are 
devised to reduce the available information to some few key numbers. 
The theoretical background of the descriptors used in this subsection was 
introduced by us in a previous work |85| . 

The shape analysis, based on the eigenvalues of the matrix with the 
moments of inertia and whose results are shown in Fig. 01 separates the 
clusters into being overall spherical, more cigar-like shaped, or more lens- 
like shaped. One can see that only few clusters have a spherical shape 
(these are found for the energetically lowest isomer for TV = 4, 6, 13, 17, 
34, 38, 54, 79, and 140, and for the next one for A^ = 42 and 116), aU 
of them corresponding to high-symmetrical isomers (cf. Table ^ and, for 
the lowest-energy isomer, most of them to the class of magic clusters. It 
is interesting that the average value follows more or less the same curve 
for all the three isomers, with some deviations at A^ = 130, 146, and 147. 
Also the largest differences show a similar behaviour, except for some few 
cases mainly for A^ below 40 and between 80 and 85. Therefore, except 
when the eigenvalues are all very similar (which occurs for A^ around 50, 
70, 100, 116, and 140), the overall shape (i.e., lens- or cigar- like) is the 
same for all three isomers. 

The construction of atomic shells can be easily seen from the distri- 
bution of radial distances (i.e., the distance for each individual atom to 
the center of mass) shown in Fig. for the ground state structures as 
function of the cluster size. Up to A^ around 50, no trends can be iden- 
tified, with an exception around A^ == 13. But for A^ just above 50 a 
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clear tendency towards shell construction can be seen for the first isomer. 
This corresponds to the formation of the AU54 icosahedral cluster. Also 
for A^ close to 110 and around 140 shell constructions for the lowest-lying 
isomer are observed. In the latter case, this corresponds to the formation 
of an octahedron. The radial distributions for the second and the third 
isomers are not shown, as they are quite similar to that for the first iso- 
mer. Particular shell constructions are found only for highly symmetrical 
clusters corresponding to A^ = 42, 48, 80, 101, 116, and around Auiao for 
the second isomer, and around A^ = 40, 60, 116, and 130 for the third 
isomer. 

In Fig. 1^ we show the average and minimal coordination numbers 
and the average bond lengths of the clusters. We define two atoms as 
being bonded if their interatomic distance is less than 3.49 A, which is 
the average value between the nearest-neighbour distance (2.89 A) and 
the next-nearest-neighbour distance (4.08 A) in bulk Au. Moreover, we 
distinguish between inner atoms with a coordination number of 12 or 
larger and surface atoms with a coordination number less than 12. 

Fig. Eta) presents the average coordination number as a function of 
N. A saturation towards the bulk limit of 12 is seen, although one has to 
remember that even for the largest cluster of our study 94 out of 150 atoms 
are characterized as surface atoms. Also, the function increases in general 
with the size of the system, with oscillations in particular for the clusters 
with N = 17 and 18, which is due to the formation of a tetrahedron for 
Aui7, and a structure with C41, symmetry at A^ = 18, respectively. The 
latter has already earlier been obtained with the EAM method (see Ref. 
|74j). but it is the first time that a tetrahedral configuration is found for 
the Aui7 cluster. 

The minimum atomic coordination for each cluster size is shown in 
Fig. |BJb). The existence of low-coordinated atoms, i.e with coordination 
numbers of 3 or 4, could point to the occurrence of a cluster growth, 
where extra atoms are added to the surface of the cluster, whereas higher 
coordination numbers could indicate a growth where atoms are inserted 
inside the cluster, or, alternatively, upon a strong rearrangement of the 
surface atoms. The latter is the case for the gold clusters, with few excep- 
tions at A^ = 14, 17, 18, 78, 83, and 134, where lower coordinations are 
found. The lowest coordination corresponding to Aui4 is in connection 
with the formation of an icosahedron plus one additional atom on the sur- 
face. At A^ = 17 and 18, some structural changes take place, as discussed 
above. Auyg and Auga correspond to structures with a decahedral motif 
capped with one additional atom. This is also the case for AU134 where 
the C2« symmetry of the decahedral structure corresponding to A^ = 133 
is lowered by the addition of an atom to the surface. 

Fig-EIc) shows the average bond length as a function of the cluster size. 
The dashed line corresponds to the bulk value of 2.89 A. The average bond 
length for all the structures is smaller than the bulk value, especially for 
Aui7 and Auis, where more compact structures are formed. However, this 
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property approaches the bulk value faster than the average coordination 
number. 

One important issue in many of the molecular dynamics studies on 
gold clusters is to identify how the clusters grow and if the cluster with 
N atoms could be derived from the one with A'^ — 1 atoms simply by 
adding one atom. In order to quantify this possibility we use the concept 
of similarity functions, introduced by us earlier |84II85| . 

The similarity function S, shown in Fig. \7]i^)i approaches 1 if the 
Avln cluster is very similar to the Auat-i cluster plus an extra atom. We 
see indeed that for N up to around 50, S is significantly different from 
1, confirming that in this range the growth is complicated. The most 
pronounced peaks occur at 6 < TV < 9, 15 < TV < 20, 34, 38, 39, 52, 
56, 79, 80, 85, 111, 126, 140, 141, and 145 < iV < 147. Many of these 
correspond to highly symmetrical clusters, however some of the clusters 
with larger peaks {N = 39, 56, 62, 85, 111, 126, 141, and 145) have 
lower symmetry. The octahedral Auss and the low-symmetrical AU39 are 
structurally very different from their A^-1-atom neighbours. Ause marks 
the end of the icosahedral shell built between AU52 and AU55, and the 
clusters resume their disordered growth. The octahedral Augi is followed 
by the disordered Au62, and the decahedral Augs comes after the disor- 
dered Au84. Between the decahedral Auno and Auii2 lies the disordered 
Auiii. The addition of one atom to the disordered AU125 leads to the 
formation of an unfinished but regular decahedron at iV = 126. The dec- 
ahedral Aui4i comes immediately after the octahedron corresponding to 
N — 140. Although AU144 has partly decahedral construction, its iV-|-l- 
atom neighbour is disordered. It seems that for each cluster size there is 
a rearrangement of the gold atoms, and no particular growth motif can 
be identified. This, in turn, means that the cluster growth is very com- 
plicated and it is difficult to consider it as an one-by-one atom addition. 

Finally, some selected, high-symmetry clusters are shown in Fig. |S1 

5 Summary and Conclusions 

We have determined the three energetically lowest isomers of gold clusters 
in the range 2 < N < 150 by using a combination of the embedded-atom 
method in the version of Voter and Chen (for the calculation of the total 
energy for a given structure), the variable metric/ quasi-Newton method 
(for the determination of the closest total-energy minimum) , and our own 
Aufbau/Abbau method (for the determination of the global total-energy 
minimum). Although the calculations provide a large amount of informa- 
tion for each individual cluster, instead of discussing each cluster sepa- 
rately, we focused on identifying general trends such as total energy per 
atom, overall symmetry and shape, average bond length and coordination 
number, and similarity with N — 1-atom clusters. 

The version of EAM used in the present calculations is parameterized 
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to bulk, as well as to the dimer properties, which allows it to describe 
properly the properties of the smaller gold clusters. 

This study predicts a number of particularly stable clusters, i.e., 'magic- 
numbered' clusters that in many cases are in agreement with results ob- 
tained by first principles and other semiempirical studies when such exist, 
but the advantage of our study is that the structures were obtained by 
using a completely unbiased approach. These magic numbers were clearly 
visible both in the 'stability function' and in the total-energy difference 
between the energetically lowest and higher- lying isomers. 

We also found that even for our largest clusters the binding energy 
per atom has still not converged to the bulk limit. Similarly, the aver- 
age coordination number is far from the bulk value, but higher than for 
nickel clusters, where several structures with shell constructions and cor- 
responding low coordination numbers were formed |85j . The average bond 
distance for gold has not reached the bulk value, due to the rearrange- 
ment of the atoms for each cluster size that leads to the formation of very 
compact structures. 

The shape analysis showed that roughly spherical clusters corresponded 
mainly to the energetically lowest isomer, but in some cases also to the 
second-lowest one, and that these often belong to particularly stable struc- 
tures. 

By analysing the distribution of radial distances as a function of the 
cluster size we could identify a region with N around 55, where a shell 
construction was formed. Comparing to previous results for nickel clusters 
|85j ■ where clear shell constructions were formed at N around 13, 55, and 
147, here the atoms rearrange for each global minimum, and therefore 
particular shell constructions can not be observed. The similarity function 
also points to the lack of regular growth. 
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Figure 1: On the left we show the binding energy per atom as a function of 
size for the energetically lowest isomers of A^ up to 150 with the dashed line 
giving the bulk value; on the right is displayed the difference between the 
total energies of the lowest-lying isomers with up to 80 atoms obtained by 
Sebetci et al. and those found here, using different structure-optimization 
methods but the same embedded-atom approach. 



19 



Table 1: Point groups of the optimized gold clusters. N.k marks the 
energetically fc-lowest isomer of the Au^v cluster. 
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Figure 2: The stability function as a function of cluster size. 
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Figure 3: The total-energy difference between the two energetically lowest 
neighbouring isomers as a function of cluster size. 
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Figure 4: Different properties related to the eigenvalues /„„ of the matrix 
with the moments of inertia. In the upper panel we show the average 
value together with points indicating whether clusters with overall spher- 
ical shape (lowest set of rows), overall cigar shape (middle set of rows), or 
overall lens shape (upper set of rows) are found for a certain size. More- 
over, in each set of rows, the lowest row corresponds to the energetically 
lowest isomer, the second one to the energetically second-lowest isomer, 
etc. In the lower panel we show the maximum difference of the eigenvalues 
for the three different isomers. 
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Figure 5: The distribution of radial distances (in A) for the fowest-lying 
isomer as a function of cluster size. Each small line represents (at least) 
one atom with that radial distance. 
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Figure 6; (a) the average coordination number, (b) the minimum coor- 
dination number, and (c) the average bond length as functions of cluster 
size. The dashed lines in (a) and (c) show the corresponding bulk values 
for gold. 
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Figure 7: The similarity function as function of cluster size. It describes 
whether the cluster with N atoms is similar to that of iV — 1 atoms plus 
an extra atom. 
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Figure 8: Some Au^r clusters with high or pecuhar symmetry. 
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